


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1989 


Inexpensive global location and tracking 
systems using geostationary satellites. 


DeFries, Danny L. 


http://ndl.handle.net/10945/27094 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist sia Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ia) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


» uw 
















pores Reet We) 

2or AIR SA tensor ated emai 6 A Qian area 8 08 Oak 

Mohd MLL peBbwot® Ded an Ocoee weir Wwe ATH err Ty 

ACTARAET DEA WRAP Wg A 1 1ke TN Or Reargancass tian D moee tRaBentoh spits ofa. outed BAAR: (4d oRetie 81011 SASS See a eee 
pars 















4 i Bae ioe tk Apapeanstsn Caek Wale eon Arar Bet Map. A Qs, and alas oAeOsa: me AcRsOGes Aghi Qotnaoer 5 ates 

rs a 2 areata. maker ara nth Jinn yan hs feast So RAR Dam As @ MARA Neg Laece—n 8s gh RATA Ogs mancnn® AT RNs = ReneS et aay 

A 5 ‘ your at Agemree 4. sm nstanea Resear OAPh: AEsK 2 REAM In AIT ERIPALRS. DAA ABOODED MgC rhs pan Raa FAO MERRE TERR Rar 
‘ can faa ane imams Be! WT BALAN Gy earn Pera dag aah MAb Ry © EERE: OARSTANE 18 1a OA REYES OEE fe Oana TT 








sara. 










ART eAtad, op hgba ante: Raat OCR meh SA RAR MEU AD, 
Mito’ gigemanrece s af Alanehaece Aotnes Rn Reasenasene AANA. 9 ge nar RL Sn QAEER) argne ee eT p eS at Ameen ee 
Mi aa ray DeRAAMINgARRTD. Roast 19: Leeery peppers 











OF AR) CA GRER SHALE TAM FEMA A MOORE IATA RED a: MR NOT TAN 








































i“ Prag pehsaapwrcyts th ored opty oar per kig phere p repen aberpear my at = earaterensonpengend 
sore ee am forest sheng aphpe pane Var -te2onae, ona a ete id ARRON pow Janennes ane am eeeE oS Se 
. ‘ » Telartheyt peer aim. Sfaeraakgadoligauspnaas Senghdnuncens;=;osanue sora avapars oe -apu/ernenee Moweedtenaa nafmeararenren-ae meme 


2) pm Daj Ahm tna aOnay mcarhyigheks RDSE RAO a dad 
fonsbiA ras bebacawasaontaynat egy nonat agantia sbanoeeh reat 

res os actus rer ant Seale nenprrtiie crit epsrprpem 
naan smopapist bsat@usnee.ssvancasane rao: igarorharan 
Avid teuen on im, Wh a saachelissacepunatireeahe W585 Aan. ates 


soa dls, Be WRAL CAR 

ROBIE AAS? gh deste Klas 
Rares, kcshate RMynce naman, Men, 
AMAL Rated HEE featce dn Ns. 
























Stl Ap tare 0 an ADAM SD e RN MME 
erased Urokne toe cere Wire 
Beinn ge UE Mey 
+ PR RO ARAR Mm Nagne AMA AREREN Ma Bum nt 





4 PRE AL Pe Ay RoW PY 25848 Sowmnse ANOS NRAATECS, SALEM AMO! Dy Rhy Ret 
te sion aes = 















eeania 
rear ny 











CARAS ES faite bs WU Une ae 5 aves ba urea 
eer errr’ 4 “ansm tne “Are aN. REAM ts UP ALstpeateaes Weg tomey 

fr yene 1 aNabNhs Me ae 9 Mananons ssen4s@aactnbor ts AAasbund ty. Sal A86 An srt AWtS af mpnamans ot 

ates ase Peter ae hate oy ARNDT ADe RRS Me Hy Tani ALNGG' Biase Mian Ae be nirag YM Hi 

Am Moms ts Maeda Ae MS me. Ae Ant thay BD DedEN OLMEAI AE ORL era agN pose eine 1 Teh aibulR O39 0B. aman | Sp PGAA ArT IRyRy Ame ae Rte 

AIBN As Dot em Ta SIA BAP) pA Late eh eltiae lnk ted fy tab yee ee ee a 


1GAGAE +e Ra tanaes . ns” ge staoe ¥ 
ET ¥ Bunghn Aas Bn SALA eed AAD Lela ee so pit intel ia el chhie arin iaetaee anid orer 






























papeeents 
bee Ree 4 Yee RJR PRS UDNOE ABO oy Uy. Mekeaan shen Lampavac@e, ha task AY Dos Benet 
. . ASRAM AMM ABANMA Mion dh) OL Pada le ois Mi 8ee 70) Se ee Nik cia Tecemens ma cadets nan SAE es et 
« . FD Mhed Sc BV Aes sR aAyIng bales OO: Aa ML SOND. DADE i es em, BR 
‘ ii 2 AAOANG ENTRAR TR Akg We Ary ER BEAN eee ayy err ery 
ed . Suarnuny at? 19090 PE ETNA TAG POSER H hom ASR, Peet: 
(<2 8 nia kad Arnal 











ha 
Ys Fires 
























2 eee Safad Leama Din ar 8 medamstinel 
Sage vee SRG ahen we dL LTA Fun ibe. pres 
. Va die oe te Ren eR y ims AStaseaetamrae ate 4 a getn aah mls A iMNGE ort Age a 
hae Vena gt B20 aboas Dota ayer teary 
ne tee 88 aas teams Ani acs hoe ASAZOI: raged 
5 ‘ 4 > ape Ge arshtia wee 





he ang ga nits 
ror 5 















ASR Eaaisare oy gona 
DS teak 






oO ae 
ec Sdrsehngcsegnesrenanta 






a 































is ' . . DASHING deen Art oey AEA peepyms tn staan tok. Mpqun ia Rear neRe ies we ee eon OR oma a a 
. oe : Jt Saecete Kingdon Tein iearanae gael wrctreta nh iafocoes ot sea aBearapader ates anaes wtbrsetae —tbmn 3 erento epee’ 
. eo ats a yt ayn a > Speer 2 ty i beast RaliPAaQye BID AB AEROS Eee pam pmamon, erage Bon mpase se awee! jerome tee wtadin 
ery oe ay re amr CL Ly peal weasel pm Leon Sang be etic i ALR Rpg, mgd van anagem tna dt Al 04 
‘ ia 

















oD ies dean Fombeline ae el pachag ata ena! whee oo 
eee midatenn seabed 3 







ek Ne! He 

a rie 

* Togeg pest aowe . 

roar] ad” Ra UNS punter tM pauls Om TaRoy * 

Vin N go agd aamiuie pega, Cuegaminyysed he Rage) ostsmmed. 
Yai 

















ae 
: . sie tgs st anaek ie oN NC ert ie nemesis 
com: in 5 © Som be hep Re 1gtehl htee eo obits ore mimer nner iy 
' ft : c &  rgStophosa “upegens Sjratae whe Late’ Poet perereseteribiy b prc etrits fem te tren be fe 





TRVeN mayy Sq UANanpOpen Simei ms Costes @ 








: 1 aemtgege mcstp Ake! 4b fo pepe facesSp tmestspiae ne gi ugehctar mmm ate eek en he 

te Feb Wane BOLE br Mase ding Rr ia 1 ty Bs ' Dna arama Se oread af eae 
rskonaasigaeilen = fase Payee Te ain yy SSS aes =a pre pb Te Reese se ae ramon ee on nae 
ee Fy Br ron Tsasans, oe Sheep yeaa sty (dang Lece<hpepence goulleh wean an 
sclavgrsAca td ieortalneandsckeeney cotraghs e aahage cave pectacqnesteadhs aft tees aeeretarer cos 

cheatin band samt Fone arent seabed nd tas ani ans SI 
Plincecneect Ts sconpeeannae ty Ms i 
Date ot ae reccese .F re 

Peary ek, oreo 
fog idager Setanta pt lay Fe 
Thee, PALRISE weap be 













































2 
Cer Sake. 
metre be 













ogee Bey ats a0 Situmes aeene tamer he! os 4 
an, CM genic er een eee ee pee en 
i. Tas A gle was epte ee 









Ui cotehe 


ip t Tigiiet a Lees. 






















re Te 
aR RATE RES 5 










pert yi 


y E t . y Ato err a Ae 





Pore 
. 


es 

ep deavreney bee 
Sra key” he ens 

se dee Sata wl 
MPEG tes lear cars Brat 

Fed acti Newer 

APPIN 8, of SEG 











































pea re eel tea 
Wa as RE ety 






mY 
na. 











‘eo ae, ay 
: aligns 4.032 
3 er ie 





































eure 
1, vee coer 
: Cee ss Beir ite 


SA se ae 
waleratbet Seth 
i 

fyobarin 
ay ae Sees 


“ ye Pap ee 
CoE ew ba Pong cyreey roel 8 Conkeapepenepres wamyrege ee ae, 
hal cae ollie Rl od abel Die ion Ports pegel area Ke 


Sp sgyetel pee oper eet 
oT apt pre fore Sree erica 
Cathe bat Lace acshes weppl mghubee caplana™ 2 x ei 
[Se 5 ucoen t ptinenetndedege se Seltoters ux Serre nae 
eg orartig on mpage Spyapenranetne are egg SSe- 
eer ee aman en yee se pene eee oe “a= ae 
ab ati Fernie mite ap epaaeey panentoer sagan Amen camer aes OS Br to 
Cn Rint che dg cna epee a hag en ne a 
Scoriemtoneen A Senn eigdt aes wim wens tien am Sens 
; : adalat pence ood ggg op gee 
pastel Sten be m-? eae eym spi= seeberes—- =, naweplet in sph don v5 Dee eee 
“ F fol Recap katate pe a sinatra gorge nas ne Bar Mesos ge a nar c aeat = om aes SRN 
: 5 3 3 a P Joy spate “oon pint ytcte ree i whe Sorin gion ene prenacaonet eben ciel gigrage eae mis a nae arene 
een poe ih oe Cig rag wa eeyuin eens ean ee ee pge ae nee nn ae 
a op fo. aL gay goto prtepiee papi paeinwis perretll ome = ares 


oe ig py ense » 


mo 
Se eerie eee 











vi agit wt 


Aerie paceere 























“. mee ee &; 
eps, ep ey es 
Saag sto eeepc sf 
peat eee Bigena ~eoeastinpe aon gage =o 
wo Ae EE Sep ecicpepme aor oe 














a: Pe: 
Bee rie brett ego 
ae enn aman poe * 


ape cege aoe OS 
Ee OS am are = neat 


ee a ee giana ope 





2 
of 












my om 
eee 






































ais pte 
Serntrrn 


Ponte 


colle eines Es : 
Meee rua orton eepeetoneiored 


ee areal Geer eal meen me 
PM [eer eoe gates saree § yoym ange mbemmeyergicuee Mew hays egmenne plemyS~e yer ymaree wee . 





















oa ae ic tape siete ap mre epee One oe ae § cash as ana 
Scan Saree ti Epa septa create he ete ee heen eee et 
3 G z mas roca eety wiste ait bat Sccern trey dahon beg eeate ce eters patel years oteeg nas peabana tiem yaa ona ab ity on acagia eet sop nae oo igeee cea 
ELS aden gene yen thee Seg stg a pk | Feb r deveracra-enpealalies tapas ahem ae 9 mgeene ewan renee pel golf! aie imran araae 
; : San | ination Rete 18 Tay ates pays aveenigemey crane ween gem erat ap tomer gegen oni eae 
: : an ; sielewcep sl pedaaeter wptetelgyl eect ten meg -yeernae porch Sp Sater ew sone proeraar 


Ree aeeny y 
efit cd cdsrbiane 
erty dae peor Seer 


Soran: 






avers yea 1 Fo ol reper sen e os 
erp tele eletincctes © fea aarepeers beaten gare ona 


Degree ease pemete yen enerereth 2.” pcamantat 





























é : ee oe rgd es ever wb ¢ - = 
: Se eras ee eee ren eh ce a 
e “ 1 Pant ete eae eterene ha = Oe are was Se eT demeinegi 
a ES Na Nepal nny nce Spleen petra et 
: . s oa fe papi ics 








ghee pes nee ctereme ae tree ee nee eel ae cir enna 
Concer setae seein sree ieee eee eet Rae ses age manana’ 
er Ye go Dg ik a tn part pos PRA on partment 
A Sa Parra 4 Po esa yteset Weryne/FMaTED “Wye aay caeny Wren) Mpeverinere Tv eee ona ae epee a we 
io. re a9 FO a Pe BELEN a em INTRON NIN ET Pe 

Le UR UPEU ey Bagt DUcvte=e= Web Ce. oNCBEETUREEH: OI NpeeySIyeweeen oo taney meee ate On uhyemn bees — coe meee meet 


pee 














are rouge ryan, Bemeiryraern” 
i 7 Menrere diet pemeurvenscr| Pp Aare opel to oe vey haprawye bee Wiyt ovate terp aevece apse ere AO 
yen hemesoner Wat Sa atop bo mee Fk ye genet MUN RED ow en eee F Rare te Ae aee ee A 





Sg? OOF Fi eg ere PR ee TTe Ls 4 

- . ° Ae sul warpadneenki nee ee patsle tore ermeuia pactoeeee 
Tras tue. -one OA Py LET Di at itcre byes whete eraser matey one y-T care 

ar) ia ee FE a ee ech heave 

tor ae posta osteeye worn 

Rt ealiadhe 










































a 
ra ja # fr jem Ta a Were 
Pabrast ante cer 8 Ut oto ruiaeeta pe artery ee eee 


TE ease y anrupgietatvanihe sar teetee emey oy pnievewt@nreereee ee eew rarer etImeyemre grmeateapearpert eee aE 
¥ saw ThE aS BEE NEY NET ETTAe SW ES” Re org ee eee ap pres See wares 
Atte ees ds “pene Se aiee 


NO pegs a8 a ENT 
Wastes aia ele trent pte Gonteayeyowes waay apres ery were ee Fat Ae Bares rie were 





© eis tee frit: 
ry ere wigs * 


7 
TH 
or 
tf 
CM ATE rs tel oh IM st ae tense 






YMA ae ct A Acad ee lot le 















































































aa eae Ps aeate> meet eT 
. . a . “y 
: . Arers 99 pEF ony ‘ 
a ; ‘ ey wean Sy Bea rat 
tae tike sewn pie Sete hea Hy tera 
_ Meter eppnitiner ce ee pinaergcusye tare -y aver fst pare ORe ba AAA 
1 L ed rT ~ *, ‘ 
: eer PpetariiocG ptate s ctovieye aretonaspst? -feeeesw sau Crmgons peta park pevvartre ey PATTaNs 2 repent 
PS De Mga cl a tale Laie oid Al ade ha da 2 Peele eere: aroun a erarerrere 

ye rycen 0 wre a yey Gee eer gE? ret vs Pei aag teeter et caces oe 
: a petri del age ae remy teen rctinimek toons orton! pineme nd Me chenronty-stnedea-aph 

F Sao paige aepage olka aoe atin 
Bess : AEs pate pte one gene eeaventee eae or SSetarenivamictentte act enw aiytvopre “aneareen ovone: wow entenen 
eat kere hat eae et ota en Rd ee Atlee asin AE Y patie vthrsrdntereo re 
* % . Pe ie et eee es tre ee re oe lathe et ad Neth satel eG dad a cchdentesethamrenacuietn 
j . . ry ” Feri Tie tyke Satan gre essere ecarere etege gen lorarash Cetra Sie a PA ame pone genre 
rf e bees leans poise 6 Fete Fa bucgematy 1) bra warmest arate erage p coRiO Kate arses he Ie weer ebrar er every y< CUR PNd! <Niwe soaves eecieewmeenre 

5 oe us + orp ygetey Lyrae Vaan: pepe ert Mam Aiewed baat speicr ttre at bey 
4 . . < Je) a3 Dhue ees = Ee aatyeet dene gotten Dwrerknegt eee e aera or 

é Pr ‘ Behn ey Mate es He Sie wieere ae 2 cmgins porta ee 
rushes yo geroiyaata Ps: ig rs ERMA OE Pn WN! 8 ana eR eae 
ay oe PP el thteoe Sir nde Hah sermons (Prine cme cemewislaneeae 








alts 
qoreemnit we 
OER ET pried ivee 














NAVAL POSTGRADUATE SCHOOL 


Monterey, California 





THESIS 


INEXPENSIVE GLOBAL 
LOCATION AND TRACKING SYSTEMS 
USING GEOSTATIONARY SATELLITES 

by 
Danny L. DeFries 


oer 


June 1989 


Thesis Advisor Tri T. Ha 





Approved for public release; distribution is unlimited. 


TOLLO6A4 





aclassified 


lurity classification of this page 





REPORT DOCUMENTATION PAGE 


1 Report Security Classification Unclassified 1b Restrictive Markings 








1 Security Classification Authority 3 Distribution Availabilty of Report 

) Declassification Downgrading Schedule Approved for public release; distribution is unlimited. 
Performing Organization Report Number(s) 5 Monitoring Organization Report Number(s) 

h Name of Performing Organization 7a Name of Monitoring Organization 

‘aval Postgraduate School (if applicable) 62 Naval Postgraduate School 


: Address (city, state, and ZIP code) 7b Address (city, state, and ZIP code) 

fonterey, CA 93943-5000 Monterey, CA 93943-5000 

h Name of Funding Sponsoring Organization {8b Office Symbol 9 Procurement Instrument Identification Number 
(if applicable) 


: adress (city, State, and ZIP code) 10 Source 110 Source of Fundmg Numbers Funding Numbers 


| 
Title (include security classification) INEXPENSIVE GLOBAL LOCATION AND TRACKING SYSTEMS USING 
iEOSTATIONARY SATELLITES 


2 Personal Author(s) Danny L. DeFnies 


3a Type of Report 13b Time Covered 14 Date of Report (year, month, day) 15 Page Count 
laster’s Thesis From To June 1989 43 


} Supplementary Notation The views expressed in this thesis are those of the author and do not reflect the official policy or po- 
tion of the Department of Defense or the U.S. Government. 


7 Cosat: Codes 18 Subject Terms (continue on reverse if necessary and identify by block number) 


ield Geostationary Satellite Navigation, position determination. 






? Abstract (continue on reverse if necessary and identify by block number) 

Inexpensive Global Location and Tracking Systems are currently being designed to provide the civilian market low-cost 
adio position determination. This paper discusses two possible designs. The first design employs 3 or 4 satellites, depending 
n whether altitude is known a prion, each transmitting continuous ranging signals. The user transceiver receives the ranging 
gnals, measures the time differentials of the receipt of the signals and transfers this information to a control station via a 
atellite link. The control station computes the user position from this data and sends the position coordinates back to the 
ser via another satellite link. In the second design, each user transceiver transmits a unique code to the control station via 
ie 3 or 4 satellite links, again depending on whether the altitude is known a priori. The control station measures the time 
ifferentials of the receipt of the signals and determines the user position. This position information is then transmitted back 
) the user via a satellite link. 


0 Distribution Availability of Abstract 21 Abstract Security Classification 

SI} unclassified unlimited D) same as report OC DTIC users Unclassified 

da Name of Responsible Individual 22b Telephone (include Area code) 22c¢ Office Symbol 
m t.Ha (408) 646-2788 62Ha 


D FORNI 1473,84 MAR 83 APR edition may be used until exhausted security classification of this page 
All other editions are obsolete 





Unclassified 


or 


Approved for public release; distribution is unlimited. 


Inexpensive Global 
Location and Tracking Systems 
Using Geostationary Satellites 


by 
Danny L. DeFries 
Major, United States Army 


B.S., University of South Dakota, 1975 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENGE IN ELECT RICAER ENGINEERING 
from the 


NAVAL POSTGRADLATE SCHOOL 
June 1989 


ABSTRACT 


Inexpensive Global Location and Tracking Systems are currently being designed to 
provide the civilian market low-cost radio position determination. This paper discusses 
two possible designs. The first design employs 3 or 4 satellites, depending on whether 
altitude is known a priori, each transmitting continuous ranging signals. The user 
transceiver receives the ranging signals, measures the time differentials of the receipt of 
the signals and transfers this information to a control station via a satellite link. The 
contro] station computes the user position from this data and sends the position coor- 
dinates back to the user via another satellite link. In the second design, each user 
transceiver transmits a unique code to the control station via the 3 or 4 satellite inks, 
again depending on whether the altitude is known a priori. The control station measures 
the time differentials of the receipt of the signals and determines the user position. This 


position information is then transmitted back to the user via a satellite link. 
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I. INTRODUCTION 


Today, the United States Department of Defense (DOD) is utilizing the NAVSTAR 
Global Positioning System (GPS) for accurate position determination. GPS employs 18 
geosvnchronous satellites within 6 different orbits with 3 satellites in each orbit. GPS 
provides highly accurate position determination to an unlimited number of users world- 
wide, although the position determination accuracy provided to the civilian community 
is intentionally degraded due to national defense reasons. The degraded mode still ena- 
bles the civilian community to attain accuracies of 100 meters. 

In addition to the high number of operational satellites that are required to support 
the 24-hour world-wide GPS system, the user equipment cost is also quite high due to 
its complexity. Even though the user equipment cost has decreased considerably since 
GPS was first introduced, the civilian community can not afford the user equipment. 
To lower the cost of the position determination system so that the civilian community 
can afford to install position determination equipment in one’s vehicle, a simpler system 
has to be designed and the total number of satellites deployed has to be reduced. 

This thesis presents two low-cost concepts of vehicle location and tracking using 

eostationary Satellites. The first concept is borrowed from the concept of geostationary 
satellite navigation systems presented in [Ref. 1]. The second concept is a new concept 
that is presented in this thesis. For the systems to be utilized by the civilian community, 
the systems have to be lower in cost than GPS, provide reasonable position accuracy, 
and be reliable. To cover the United States, each system requires four satellites to de- 
termine longitude, latitude, and altitude, or three satellites if altitude is not measured. 
Like the geostationary satellite navigation system, the location and tracking systems will 
also deploy three or four geostationary satellites parked in orbits on the equator at a 
latitude of zero degrees and altitude of approximately 35,786 kilometers. The difference 
between the proposed location and tracking system and the geostationary navigation 
system is that the location and tracking system user does not determine his position but 
relies on the central control station to determine it. The satellites continuously and 
synchronously transmit a satellite-unique ranging code to the user on earth. The user 
receives the satellites signals at different times due to the path length differences between 
the user and each of the satellites. The transceiver measures the time differentials be- 


tween a designated signal and the remaining signals and transmits this information to a 


control station via a satellite communications link. The control station calculates the 
user position from the time differentials and transmits the user coordinates to the user 
via another satellite channel. By placing the computational burden on the central con- 
trol station, the location and tracking transceiver can be made at a lower cost. 

Position determination in a geostationary coordinate system is normally defined in 
a geocentric coordinate system as shown in Figure 1. The x-axis is the intersection of 
the equator plane (0° latitude) and the Greenwich meridian plane (0° longitude) and is 
oriented from the center of the earth, the z-axis is the polar axis oriented from south to 
north, and the y-axis completes the right-handed Cartesian coordinate system. We des- 
ignate the unknown user position by (x,y,z), the known positions of n= 3 or 4 
geostationary satellites by (x,y, 2z,) Where i= 1, 2,...,2, and the known position of the 


central control station by (x9, 3, Zo). 





0° Longitude 


Figure 1. Geocentric Coordinate System:  [Ref. 1] 


The main function of the transceiver is to measure the differential distance 
a; = d, — a, i=], 2, sey nt (1) 


representing the differences in the distance to the user from one satellite, arbitrarily 
designated as satellite one, and from each of the other satellites. The variable a, repres- 
ents the difference in the distance of satellite one to the user and that of each of the other 
satellites to the user. The variable d, represents the distance from each satellite to the 
user. The differential distance determination concept is discussed later. The squared 


distance from the user position to each geostationary satellite can be expressed by 
2 2 2 2 : 
d; = (= x7) = (yy) 42 (B= = u t= ae seeg (2) 


where (x, y, 2) represent the users coordinate position in the Cartesian coordinate system 
and (x, },. Z,) represents each of the satellites’ positions. 
For n= 4, (1) and (2) combine to form a set of four independent equations and four 


unknowns, namely, x, v. Z, and d,. These four equations are 
2 2 2 2 : z 
Cee ee ey) t(2—2), ta 1,2,.5,4 (3) 


which can be solved using the measured data a, and the known satellite ephemeris 
ees tk a8 4. The value a, is found through the use of a correlation register and 
ranging codes which will be discussed later. A detailed discussion of the concept can be 
found in [Ref. 1]. 


For the three satellite system, n = 3, we use 
x ty +27 =(R+AY (4) 


Where R = 6378 km is the Earth’s mean equatorial radius, and h is the user altitude 
{assumed to be known), as the fourth equation in addition to the three equations 
(goa — (x —x) (yy) t(z—z) i=1, 2,3. 

The first system requires each satellite to broadcast a unique pseudo-noise (PN) 
code in the frequency band 2483.5 - 2500 MHz. The PN codes are selected from a set 
of Gold codes which allow minimum interference between signals from various satellites. 
The satellites employ extremely stable on-board atomic clocks which are mutually syn- 
chronized. A central control station which can view all satellites is responsible for the 
clock correction via a communications link in the frequency band 6526.5 - 6541.5 MHz. 


The user transceiver measures the time differentials between one satellite arbitrarily 


designated as satellite one, and each of the other satellites. The accuracy of the meas- 
urements and their effect on the user position determination accuracy will be discussed 
later. The user transceiver then sends the measured time differentials to the central 
control station for computing the user coordinates via a code division multiple access 
(CDMA) uplink in the frequency band 1610 - 1626.5 MHz. The central control station 
calculates the user coordinates using the received time differentials and the predicted 
satellite ephemerides. It then relays the coordinates to the user via a time division mul- 
tiplexed (TDM) downlink in the frequency band 2483.5 - 2500 MHz. 

The second svstem requires the user transceiver to broadcast a PN code to the cen- 
tral control station via a code division multiple access uplink. The central control sta- 
tion calculates the user differential distances to the satellites with the available satellite 
ephemeris and then the user coordinates which are relayed to the user via the TDM 


downlink. 


Il. SYSTEM | 


A. SYSTEM I DESCRIPTION 
The frequency assignment and bandwidth allocated for each of the satellite links 
listed below are assigned by the Federal Communications Commission as listed in [Ref. 
2] 
¢ Satellite-to-User Ranging Link 


Frequency 2483.5 - 2493.73 MHz 
Bandwidth 10.23 MHz 


e User-to-Control Station 
User-to-Satellite Uplink 
Frequency 1610 - 1626.5 MHz 
Bandwidth 16.5 MHz 
Satellite-to-Control Station Downlink 
Prequeney —5150'= 5160.5 Wirz 
Bandwidth 16.5 MHz 


¢ Control Station-to-User 
Control Station-to-Satellite Uplink 
Frequency 6524.5 - 6530.77 MHz 
Bandwidth 6.27 MHz 
Satellite-to-User Downlink 
Frequency 2493.73 - 2500 MHz 
Bandwidth 6.27 MHz 


¢ Control Station-to-Satellite Command Link 

Frequency 6526.5 - 6541.5 MHz 

Bandwidth 15 MHz 
B. SATELLITE-TO-USER DOWNLINK 

1. Downlink Frequency and Bandwidth 

The communications downlink from each of the satellites to the user occupies 
a bandwidth of 16.5 MHz from 2483.5 to 2500 MHz. Within this bandwidth, the satel- 
lites must transmit a pseudo-noise ranging code needed for user position determination 
and provide the downlink for the control station-to-user data link. In System I design, 
the satellites must employ two satellite-to-user channels, one for the sole purpose of 
transmitting the unique pseudo-noise code and one to perform as a relay for the data 
link from the control station to the user. The pseudo-noise code will be transmitted 
between the frequencies of 2483.5 and 2493.73 MHz in a bandwidth of 10.23 MHz. The 
contro] station-to-user data link will be transmitted between the frequencies of 2493.73 
and 2500 MHz in a bandwidth of 6.27 MHz. 


2. Pseudo-Noise (PN) Ranging Code 

The PN ranging code transmitted by the satellites is received by the user and is 
used to determine time differentials. The accuracy of the measurements of the time dif- 
ferentials used to determine the user position is of prime importance to the overall po- 
sition accuracy. The more precise the time measurement is, the more accurate the 
position determination will be. GPS uses code chip rates of 10.23 Mbit/s in its trans- 
missions to the user. For our example, a given code length of 2! — | = 524,287 chips 
and a chip rate of 10.23 Mbit/s results in a code period of 524,287/10.23 x 10° = 51.25 
ms. A lower code length of 2!§ —1 = 65,535 and a lower chip rate of 1.023 Mbit's has a 
period of 65,535. 1.023 x 10° = 64.06 ms and a larger error in the differential distance 
measurement. The code tracking loop that would be used is able to maintain code 
alignment to approximately 0.01 chip. This is a conservative figure, for code tracking 
loops can be accurate to 0.001 chip. For a 0.01 alignment factor and a 10.23 Mbit’s chip 
rate, a ranging error of 0.29 meter would exist, where a 1.023 Mbit,s chip rate would 
have a ranging error of 2.9 meters [Ref. 1]. To better understand the fundamentals of 
PN ranging codes that determine position accuracy, Figure 2 shows the concept of dif- 
ferential distance measurement. To determine the differential distance, the user equip- 
ment must be able to perform a correlation operation among the satellites’ codes and 
internally stored copies of those codes. The receiver clock must be able to count the 
number of chips between the correlation peak of the first received code and the subse- 
quent codes as shown in Figure 2. Since a, may not be an integer multiple number of 
chips, the receiver may use the correlation value relative to the peak correlation voltage 
to determine a fraction of a chip in the measurement. As was stated earlier, a fraction 


of a chip of 0.001 can be measured. 


emi chip ' | chip -1 chip!’ 1 chip t 
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Figure 2. Differential Distance Measurement: ([Ref. 1] 


In [Ref. I], the PN ranging code and its effects on position determination accuracy is 


discussed in further detail. 


C. USER-TO-CONTROL STATION LINK 
1. Uplink Bandwidth 

The frequency bandwidth allocated to the user-to-control station data link is 
1610-1626.5 MHz. The user direct sequence spreads the data with a Gold code of period 
2'© - | = 1023 chips which is a standard period for navigation systems. Thus for a 
bandwidth of 16.5 Mliz from 1610-1626.5 MHz, the typical user data rate is approxi- 
mately 16,129 bits’/s when spread at 1023 chips per bit, if Quaternary Phase Shift Keying 
(QPSk) is used as a modulation scheme. ‘That is, the chip rate is (16,129)(1,023) = 
16.5 x 10° chips. With QPSk, there are 2 bits per symbol and hence the symbol rate is 
§.5 x 10° symbols a second. The bandwidth required (double sideband for QPSK) is 
twice the symbol rate or equal to 16.5 MI1z which is compatible with the available 16.5 
MHz bandwidth. The users must share this available bandwidth with each other to 
transmit their [Ds, time differentials, and other information to the control station. 

lie user-to-satellite uplink is the most critical link in the system. The user is 
not able to employ or utilize a high gain dish antenna like those used in the fixed 
stations. The maximum output RF power is also limited by the power available in a 


mobile vehicle. To calculate the Effective Isotropic Radiated Power (EIRP) ofa vehicle, 


Eq. (5) is used, where P, is the transmitter power out and G, 1s the effective gain of the 


antenna. 
EIRP = PG, (5) 


In our design, a transmitter power out of 80 watts which is equal to 19 dBW and an 
omni directional antenna with a gain of 4 dB will yield an EIRP of 23 dBW. These design 
values are selected with some judgment and are also being used in current navigation 
systems. 

The Free Space Path Loss (FSPL) is calculated using Eq. (6). 





FSPL(dB) = 10 1og( ait ) (6) 
where fis the carrier frequency, d is the slant range, and c is the speed of light. At a slant 
range of 35,786 km and a carrier frequency of 1618 MHz, a FSPL of approximately 190 
dB is obtained. 

A standard satellite G'T of +3 dB K has been chosen for this design model and 
Boltzmann’s constant k is given as -228.6 dBW;K-Hz. To calculate the uplink carrier- 


to-noise density ratio, Eq. (7) is used. 
< (dB) = EIRP(dBW) — FSPL(dB) + = (db! K) — k(aBW! K-Hz) (7) 


The EIRP of 23 dBW, FSPL of 190 dB, satellite G’T of +3 dB/K and Boltzmann's 
constant of -228.6 dBW K-Hz are entered into the Eq. (7), obtaining an uplink carrier- 
to-noise density ratio of 64.6 dB-Hz. This density ratio is sufficient to support the data 
rate that the system requires for successful operation. 

For the downlink, the same calculations are performed to obtain the downlink 
carrier-to-noise density ratio. A FSPL of 198 dB is calculated for the carrier frequency 
of 5150 MHz. A satellite EIRP of 36 dBW and a control station G/T of 20 dB’K are 
used to calculate the downlink carrier-to-noise density ratio. These values are standard 
satellite EIRP and station G’T values. 


The total carrier-to-noise density ratio is calculated using Eq. (8) found in [Ref. 


($)-[S0(8Y.] ° 


3]. 


To solve for the total carrier-to-noise density ratio, the values of the carrier-to-noise 
density ratios must be converted to real numbers, substituted into Eq. (8), the equation 
solved and the answer converted to the logarithm form of dB-Hz. 

A margin of 2 dB is further engineered into the system design to allow for in- 
terferences that would degrade the operation of the system. The data rate is calculated 
by subtracting the 2 dB margin from the total carrier-to-noise ratio, converting that dB 
value to a factor and then dividing that factor by the equivalence of the required 
signal-to-noise ratio which is set to two representative values of 15 and 20 dB for this 
system design. These required signal-to-noise ratios are typical satellite navigation 
signal-to-noise ratios; each will be used for design calculations. 

To calculate the data rate for a lower required signal-to-noise ratio of 15 dB, the 
2 dB margin is subtracted from the total carrier-to-noise ratio of 64.57 dB-Hz, giving 
62.57 dB-Hz. This ratio is converted to a value of 1.807 x 10° and is divided by 31.26 
(value of 15 dB) which in turn vields the data rate of approximately 57 kbps. When the 
required signal-to-noise ratio is increased to 20 dB (factor of 100), the data rate is low- 
ered to 18 kbps. The link analvsis is tabulated in Table 1. 

2. User-to-Control Station Data 

After receiving 3 or 4 satellites’ ranging signals and measuring time differentials, 
the user must transmit the information to the control station for calculation. The 
transmission is Code Division Multiple Access (CDMA) as required by the Federal 
Communications Commission (FCC) and is designed to provide the system with a bit 
error probability of 10° to 10°. (A bit error probability of 10-° and 10-° provides the 
reliability needed for a successful system.) The data stream must include a synchronous 
code word in the user transmission for the control station to synchronize to the trans- 
mission. Once the the control station has synchronized itself to the data transmission, 
the control station can identify the information that the user is sending to the control 
station from within the data frame. A Barker code word of 13 bits would provide the 
synchronization that the control station equipment requires and enable the control sta- 
tion to lock on to the transmission and recover the user transmitted data [Ref. 4]. 

Other information that is sent includes the user ID, the ID of the satellite that 
was received first, the second and third satellites IDs and the time differentials AJ in 
symbols between receipt of the signals of the satellites. These time differentials are de- 
termined by the correlator. The minimum distance that a satellite may be from the user 
is 35,786 km which corresponds to a minimum propagation delay of 119.287 ms. The 


maximum distance is 41,678.82 km which corresponds to the maximum propagation 


Table 1. SYSTEM I USER-TO-CONTROL STATION LINK ANALYSIS 


Uplink 
Carrier-to- Noise 
Density Ratio 


Downlink 
Carrier-to- Noise 86.6 dB-Fiz 
Density Ratio 


Total 
Carrier-to-Noise 
Density Rave 


If required ~ i 20 dB 


Then Data ae can be 57 kbps 18 kbps 







































64.57 dB-Hz 


delay of 138.929 ms. A maximum time differential between the minimum and maximum 
time delays is 19.642 ms. We require an accuracy of 0.1 usec. This allows a quantization 
interval of 0.2 usec. The number of 0.2 usec intervals in 19,642 psec is 98,210. To rep- 
resent this value, 2'* = 262,144 will easily suffice to represent the maximum time differ- 
ential. Thus there are 18 bits required to represent the maximum differential time 
measurement between the first satellite and each of the other satellites and 2 additional 
bits to identify the satellite for a total of 20 bits used. In the three-satellite design, the 
altitude of the user must also be transmitted to the control station for position determi- 


nation. An altitude of 100,000 meters with 0.1 meter accuracy can be represented with 


20 bits. For the four-satellite case, the field for altitude is used for satellite 4 1D and 
differential time measurement. In all cases, 20 bits for additional information are de- 
signed into the frame for future uses such as emergency reporting of position location 


when the user may be in need of it. See Figure 2 for an example of a data frame. 


or GENERAL 


SAT 4 INFORMATION 
1D & DT 





Figure 3. Data Frame (User to Control Station) 


To increase the number of channels that could exist on the system, Forward 
Error Correction (FEC) is used. The use of FEC lowers the required signal-to-noise 
ratio and thus allows more CDMA channels to be used. A convolution code with 
Viterbi decoding and soft decision will increase the total number of channels [Ref. 5}. 
There are various combinations of convolution code rates and constraint lengths that 
can be examined to provide the optimum data rate for the svstein. 

To calculate the number K+1 of CDMA channels that can successfully exist 
within the bandwidth at a bit error probability of 10° and 10-’, the following results 
from Refs. 5 and 6 were used. To calculate the number K+ 1 of CDMA channels that 


can be succesfully utilized simultaneously, Eq. (9) is used. 


= { ito No ¢ 
‘ — an( Nala 0) 


E 
The chips per bit (N), the design bit energy-to-noise ratio =.) and the bit energy- 


i 


to-noise plus user interference density ratio FT, J must be known to solve for K. For 
c-) 


reliability, bit energy-to-noise ratios of 15 and 20 dB were selected for the system. The 
chips per bit (N) is 1023. The value for the bit energy-to-noise plus user interference 


density ratio is a function of other parameters as follows. 
Ey 
p=0 i (10) 
No 


where Q(x) is defined as 


OG) = en eae (11) 


Jn 
x 


The code rate (R) and the channel transition probability (p) must be known to calculate 
the energy-to-noise plus user interference density ratio. A code rate (R) of 1,2 to 7,8 
was used in the calculations to determine the bit rate that will support the maximum 
number of CDMA channels. The probability (p) is a function of parameters that are 
determined by the code rate and system design. To calculate p, Eq. 12 is used [Ref. 6]. 


1 a Afee 
P» coded = i Bay, eg IF: (12) 


For System | design, coded bit error probabilities (P,..4a) Of 10-° and 10-” are 
used to insure a reliable svstem for the user. The number k of encoder input bits and 


the number n of encoder output bits establish the code rate (R), where R a The 


variable B,,,, used to calculate the channel transition probability (p) is a calculated value 
that is determined by the code rate, constraint length, and the free distance of the code. 
The variable d,,, is the free distance of the convolutional code and is determined by the 
and d, 


‘free 


code rate and the constraint length. The values Coy 


culated in [Refs. 5,6] at a constraint length of 7 for the code rates of 1'2 to 7/8. 


have been previously cal- 


To calculate the number of simultaneous CDMA channels that the system can 
support, the channel transition probability (p) must be found with Eq. (12). The number 


k of encoder input bits, B, , d,.. and the coded bit error probability (P,,..zq) are entered 


free 


into Eq. (12) to solve for the channel transition probability (p). Then the bit energy- 


to-noise plus user interference density ratio ae is substituted into Eq. (10) along 
o 


with the given chips per bit (N) of 1023, and the design bit energy-to-noise density 


ratio = of 15 and 20 dB. In Tables 2 and 4, the number K + | of channels that can 
successfully be used with hard decision Viterbi coding at a bit error probability of 10-° 
and 10~" respectively are shown. For example, given the code rate (R) of 1,2, the en- 
coded = bit k 1s" 1, ye is 36, and d,,, is 10. For a Py ..¢q Of 10-6 a channel transition 
probability (p) of .00703 is calculated. With a p of .00703 and a code rate (R) of 1/2, a 
bit energy-to-noise plus user interference density ratio <a ) of 5.87 is calculated. This 
value, the chips per bit (N) of 1023, and the design energv-to-noise density ratio of 20 
dB (100) vields a value K of 246. The number of channels K+1 is then 247. An in- 
spection of the tables reveals that a code rate of 1/2 yields the largest number of simul- 
taneous channels that can be used successfully. When soft decision Viterbi coding is 
used, a decrease in the bit energy-to-noise plus user interference density ratio of 2 dB is 
realized. Thus soft decision Viterbi coding further increases the number of CDMA 
channels that can be used which is shown in Tables 3 and 5 for 10-* and 10~° respectively. 
The different bit error probabilities of 10-° and 10-7 and soft and hard Viterbi decisions 
are used in the calculations to provide the reader an indication of four possible system 
capabilities. 

The major limitation in using lower code rates of 1'2 is that the bandwidth of 
the svstem that is available to the user is severely reduced. Although the number of 
CDMA simultaneous channels is greatly increased, the bandwidth of the system is de- 
creased. Further research is required in the area of the coding techniques to find the 
maximum number of CDMA channels that can be used with the highest coding rate to 
allow for the greatest user bandwidth. 

3. Satellite-to-Control Station Downlink 

The frequency bandwidth allocated for the satellite-to-control station downlink 

is 16.5 MHz from 5150 to 5165.5 MHz. The satellite serves as a relay for the user-to- 


control station data link. 


D. POSITION DETERMINATION 

After receipt of the time differentials, the central control station computes the co- 
ordinates of the user from the measured time differentials, and the known location of 
each of the satellites. 

The x, ¥, Z position determination of the user is defined in a geocentric coordinate 
system as shown in Figure 1. The x axis is the intersection of the equatorial plane and 


the Greenwich Meridian plane and is oriented from the center of the earth. The y axis 
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Table 2. NUMBER OF SIMULTANEOUS CDMA CHANNELS WITH FEC 


(Hard Decision) at Coded Bit Error probability of 10-° 


Code 
Rae 1/2 2j3 3/4 4'5 5/6 6/7 18 






CDMA Channels K 
E 
for aa = 100 or 20 dB 


CDMA Channels K 
iE. 
7 = 32 or 15 dB 


Sis 






for 


completes the cartesian coordinate system with the x axis. The z axis is the north to 
south axis connecting the two poles. 

Although actual distance is not measured by the user equipment, the distance can 
be calculated from the time it took for the other satellites signals to reach the user versus 
the first satellite signal that was received. The time differentials are converted to differ- 
ential distances by multiplving the time differentials by the speed of light (c); approxi- 
mately 3 x 10 m’s. The accuracy of the synchronization of the satellites transmissions 
plays an important part in the users position determination. If the satellites do not 
transmit at the exact same time, the time differentials will not be correct and thus neither 
will the distance differentials. 

The distance between the unknown user position (x, y, Z) and each of the satellites 
(x, 42) wherey = 2,4... can Dercaicularedeinom 


=(x—x) +p) $(2-z)? f=1, 2.640 (13) 
We further define the distance measurements 


a=d,—d, i=1,2,..47 (14) 


Table 3. NUMBER OF SIMULTANEOUS CDMA CHANNELS WITH FEC 


(Soft Decision) at Coded Bit Error probability of 10-° 


Code 
R 
_ kk eee Eee Te ee ee 





CDMA Channels K 


+ = 100 or 20 dB 


o 









for 






CDMA Channels K 


Ey - 
= S2 Or io dB 


c) 











for 


as the difference in distance from one satellite designated as satellite one to the user and 
the other satellites distance to the user. 

For a four satellite case, n = 4, equations (13) and (14) are combined to form four 
independent equations and four unknowns, x, y, z, and d. From the known satellite 
positions and the time differences which have been measured by the user, the four un- 
knowns can be solved. Sign ambiguities can be corrected as long as the user knows what 
quadrant he is in. 

The coordinates of the satellites (x, y,, z,.) can be calculated from its longitude 


eastward of zero degrees by 


xj =r cos 6; (15) 
z,=0 (17) 


The term z, will remain zero as the satellites are in an equatorial orbit with a latitude of 


zero degrees. 
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Table 4. NUMBER OF SIMULTANEOUS CDMA CHANNELS WITH FEC 


(Hard Decision) at Coded Bit Error probability of 107 


Code 
Bu 1/2 2/3 3/4 4/5 Oe 6/7 


ee 


CDMA Channels k 
for <= 100 or 20 dB 215 | 200 159 


CDMA Channels kK 


E, - 
2) rteg 32 or [5-dB 


° 





The four non-linear independent equations are 
(d,— a) =(x—x) + Y-y +2? f=1,2,3,4 (18) 


Equation (18) can be reduced to 3 linear independent equations bv substituting equation 
(14) into equation (13) and subtracting the (/= 1) equation from the other 3 equations. 


The remaining three equations are 


2 
a; ; 
= ee ery a o i=2, 3, 4. (19) 


The 3 equations can then be entered into matrix form to solve x, y, and d, as shown: 


my A2 eee) ce qs 
1 2 

Y | oy ee ees a3 (20) 
2 
qd, oe ee ee ee 


z can be calculated by substituting equation (20) into equation (18) 


Table 3. NUMBER OF SIMULTANEOUS CDMA CHANNELS WITH FEC 


(Soft Decision) at Coded Bit Error Probabilitv of 10-7 


Code 
ae aes 2/3 3/4 AES 5/6 6/7 7/8 






CDMA Channels k 
: 
for ~ = 100 or 20 dB 


CDMA Channels kK 


Ey 2 
for = 32 or 15 dB 










From equation (20) and (21), the user position (x, v, z) can be solved to obtain an initial 
approximation. The actual position must be solved iteratively or by utilizing a Kalman 


filter. From (x, y, z), the users altitude h, longitude and latitude can be derived. 





=(R+h) cos 6, cos 6; (22) 
=(R+A) sin 6, cos 6, (23) 
z=(R+A) sin 6, (24) 
where 6, 6, and h are solved by 
@,=tan“(-+) (25) 
g,= tan "(2 ot (26) 
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— Z — 
A rere (27) 





For a three satellite system, the position solution algorithm is much the same. The 
following equation substitutes for the fourth satellite equation 


x ty ee Se (28) 


As in the four satellite case, these four equations must be solved iteratively to obtain 
the users position. In the three satellite case, the user must report his altitude because 
only the longitude and latitude of his position can be determined from three satellites. 


The four non-linear equations are 
(=a) =@—x) yy) +e = ee (29) 
vty +27 =(R+AY. (30) 


As in the four satellite case, the equation with index = 1 is subtracted from the 


other two equations. 


dy = a,x + Byy + Yin (31) 
ad, =a,x + By + 2 G2) 
where 

meee ee) Le 

mere io SS a (33) 
2 a 2 2 2 
B aoa ew eo ar oe Oe ee a (34) 
= a3 1 2a) 

Pei= foeen se 
a 2a; , id 


When we solve for y im terms ol x wezet 
Y= MX + M2 (36) 


where 


Yi 8 Wy = Pa 


7 oyey? Eee > 
The substitution of (36) into (31) produces 
d, =u3x + bg (38) 
where 
ye 4. (39) 
The term z? can be expressed as 
Yad = = ae eee ae Uni x? - ite (40) 
Substitution of (36), (38), and (40) into (28) produces 
ax? + bx+c=0 (41) 
Where 
a= 1; (42) 
b= 23M, + 2y\u, + 2x, (43) 
c= uy + 2yyuy — x} — yp — (R+A). (44) 


The value of x, v, and z can be calculated from (36), (40), and (41) subject to the 
limitation x < 6378 km provided the user knows what quadrant of the earth he 1s lo- 


cated in. 


E. POSITION ERROR 

Even though the satellite 1s geostationary, the satellite does move about its nominal 
position. The drift of the satellite can be calculated to less than 20 meters from its 
nominal position. The error that results in the user position determination can be cal- 


culated. The range error can be expressed as uSer position error (Ax, Ay, Az.) 





fal, 2 3, 4 (45) 








Equation (45) can be transformed into matrix form as Av = u and with equation (14), 














x— Xy Jaa zZ— 2 
dy qd, d, 
aa) ae me 
ad,—a a im, 
yi 1 2 1 1 (46) 
2 2G v= s3 Z— 23 








a, — a, Gy = ee ad, — a, 


x= xy vrs 2 Fi 





d, — ay ad, — a% di — a% 


Ad, 
Ax 
Ad, 
eae a : (47) 
Ad; 
Az 
Ad, 


The user position error can then be simulated by letting the range error be 


Ad, = — x ei ae = 
= [@-x) + v-y) + &@ dy a—al, 2, 3,04 (48) 
ca f(x — x) =P (y—y) ate (2—z) ]2 


where 
(o% Vp z’) = (x; + Ax;, Si + AY Zj + Az;). (49) 


The variable Ad, is the difference between the distance from the user to the nominal 
satellite position and the distance from the user to the satellite position that the satellite 
may have drifted to. The variables (x’,, )’,, z’;) are the coordinates of the nominal satellite 
position with the drift distances ( Ax,, Ay,, Az, ) added to the nominal satellite position 
coordinates. The error vector v is the linear least-squares solution of Av=wu where 
v=(A7A)'A7u and the mean squared error is e=(Av—u)’(Av—u). The values 
(Ax, Ay, Az) are calculated using Eq. (46) and (47). The error contribution of the z co- 
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ordinate is much less than the x or y coordinate as the user position approaches the 
equator. This ill-conditioned problem results in a poor estimation of Az which can be 
minimized by modifving the least-squares method of solution. By first calculating the 
least-squares estimate of Ax and Ay and ignoring the Az component of the total error 


and calculating Az after Ax and Ay are found, a better estimate of the total error can be 


found. 
Ax 
y= | = (A; Ay Ay & (50) 
Ay 
bz = (Ay Ay) (Aj u — Az A, ») (51) 
where 
* a8) Bho Vy 
da, da, 
Vo xX Je = 
d, — a d, —a 
hie 1 2 1 2 (52) 
XT X3 ae ag 


Pi X4 J — Va 
d; p= a, a, = ay 
Zo Z|} 
d, 
A= Ad, 
aj =e Ad 
re en | ae (53) 
Bo Ea Ad, 
d, — @ Ad, 
La 24 


For a three satellite system, the range errors Ad i=1, 2, 3 are calculated from 


equation (45) and the user altitude error AA is calculated from equation (30) as 


7a]| 


m x(Ax) + y(Ay) + 2(Az) 


[e apr sig aP AlP (54) 
_ x(Ax) + y(Ay) + 2(A2) 
= R+h 


The user position error vector v = [Ax, Ay, Az]? is the least-squares solution of 
By = wwhich is v = (B7 B)" B’ w, where 











* — fee 27 41 
a, a, a 
XS e ge Ad, 
d—-a d,—-—a d,-—a Ad. 
p= 1 2 1 2 1 2 eh: 2 (55) 
AS ee ea mm Ad; 
x of Zz 
R+h R+h R+h 


For the three satellite case, the least-squares method is also modified to provide a better 


estimate of y. The variable vy is calculated from v = [v,, Az]? where 








Ax 
les | = Ga Bi w (56) 
Ay 
Az = (BS Be (Bw a Ba) (57) 
and 
xX ¥ = yy Ze 2 
d, d, d, 
2s = xX v = V2 z= 2) 
ad, —a a,-—a d,—a 
De 1 2 1 D he 1 2 (58) 
ames Vm emacs 
d, — a3 d, a a3 da = ay 
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Equations (50), (51), (56), and (57) calculate the upper limit of the error in the user 
position provided the satellites position errors are known in either the three of four sat- 
ellite case. The average error for various users positions has been compiled in [Ref. 1]. 
The average error in the users position increases as the user approaches the equator. 
The average error ranges from 14.2 to 61.5 meters depending on how close to the 
equator the user is. The closer the user is to the equator, the worse the error which in- 


dicates a Weakness in the system. 


F. CONTROL STATION-TO-USER DATA LINK 
1. Uplink Bandwidth 

The frequency bandwidth allocated for the control station-to-satellite uplink is 
6524.5 to 6530.77 MHz. The 6.27 MHz bandwidth coincides with the available band- 
width of the satellite-to-user downlink part of the relay. 

After the control station has completed the calculations and the user position 
is determined, the position coordinates have to be sent to the user. A TDM channel is 
used to transport the user information to the user. To calculate the total carrier-to-noise 
density ratio of the control station-to-user link. Eqs. (5) - (8) are used as was discussed 
in the user-to-control link. 

The uplink frequency is 6533 MHz with an antenna gain (G,) of 44.5 dB. which 
is a standard antenna gain for a ground station antenna. The antenna gain is much 
higher than for a mobile vehicle as one is able to use a large dish. The power out (P,) 
is set at only 1 watt per channnel since there would be more than one channel simul- 
taneously transmitted by the contro] station. For example, if 200 CDMA channels were 
used, then 200 watts of power would be required for transmission of the channels. The 
lower power out is standard for a ground station. For this link design, the satellite G'T 
is -1.5 dB K. The power out and the satellite G’T values are design parameters deter- 
mined by judgement and current systems. For the control station-to-user link, the total 
carrier-to-noise density ratio is 66.25 dB-Hz is calculated by Eq. (8). A carrier-to-noise 
density ratio of 66.25 dB-Hz is sufficient for a successful link. A 5 dB margin is engi- 
neered into the system to allow for the cross-interference degradation of the system by 
the multiple carriers that are transmitted by the control station. After subtracting the 
5 dB margin from the total carrier-to-noise density ratio, converting to a factor and di- 
viding the result by the required noise ratio factor, the data rate for the bit error proba- 
bilities of 10-° and 10-’ are 42.2 kbps and 13.3 kbps, respectively. The calculations for 


the control station-to-user link utilize Eqs. (5) - (8), which have been discussed earlier. 
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These are the same equations as Table | with only the input values differing. The con- 


trol station-to-user link analvsis is tabulated in Table 6. 


Table 6. SYSTEM I CONTROL STATION-TO-USER LINK ANALYSIS 


Antenna Gain 
Free Space Path Loss 


Uplink 
Carrier-to- Noise 
Density Ratio 


Downlink 
Carrier-to-Noise 
Density Ratio 


Total 
Carrier-to- Noise 
Density Ratio 


Margin 


; E, 
If Required "ie ag 15 dB 20 dB 


Then Data Rate can be 42.2 kbps 13.3 kbps 




















72.6 dB-Hz 











67.4 dB-Hz 












66.25 dB-Hz 





2. Control Station-to-User Data 
Using the available data bit rate, the user coordinates, identification tag and 
synchronization bits must be transmitted to each user that has requested a position fix. 
There are 58 bits required to represent the latitude, longitude, and altitude. Seven bits 
are allowed for degrees (where 2) =.128 which will suffice to represent the users position 
if the user supplies to his machine the quadrant that he is in). There are 6 bits for min- 


utes, and 6 bits for seconds, where 2° = 64 will suffice to represent 60 values for both 
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latitude and longitude. Another 20 bits are required to represent the user altitude. In 
Figure 3, a 40 bit code word is used for synchronization of the TDM frame. A 40 bit 
code word is a design parameter which is used in other communication systems. Each 
user information frame consists of two 8 bit flags, position information of 58 bits, gen- 
eral information of 20 bits and a user ID of 10 bits. The 8 bit flags provide the start and 
end of the user data. 


SYNC USER 1 USER 2 jn [=| | USER N 





104 BITS 


Figure 4. TDM Frame (Control Station to Users) 


3. Satellite-to-User Downlink 
The data link from the control station to the user is completed with the 
satellite-to-user downlink at a frequency bandwidth of 6.27 MHz between 2493.73 MIIz 
and 2500 MHz. The satellite provides relay service only and does not process the data 
that is passed through. 
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Il. SYSTEM II 


A. SYSTEM II DESCRIPTION 

System II will emplov 3 or 4 satellites in the same arrangement and positions as 
System I. System II differs from System I in the technique by which the differential 
distances are obtained. System II users will transmit a user unique pseudorandom 
ranging code to the control station via a satellite relay. The time delays between the ar- 
rivals of the codes from each of the satellites will determine the differential distances of 
the paths. Each time the user transmits a code, the code will travel to the satellites at 
the speed of light. The satellite that is closest will receive the code first and relay the 
code to the central control station. The closest satellite will provide the shortest time 
delay. The satellite that is furthest away from the user and the central control station 
will have the longest delay. It is these differences in the times of arrivals that determine 
the distance differentials as we did in Svstem I. In System I], a total user-to-satellite and 
satellite-to-control station differential distance is calculated from the total time differen- 
tial. The central control station assumes the task of measuring the differences in the 
round trip distance between itself and the user via one satellite, arbitrarily designated as 


satellite one, and each of the other satellites, We define these differential distances as 
b=(,+d))-(+4), i=1,2,..,27 (59) 
where / is the distance from the central control station to the satellite 
a 2 2 2 : 
b=Go—-x) + Go- yy “lzo a) a 2 (60) 
and d is the distance from the user to the satellite 
= (@—a) Hip ee ee lee (61) 


For n =4, Eqs. (9) - (11) combine to form a set of four independent equations and 


four unknowns, namely, x, y, Zz, and ad, These four equations are 
(4, —- bhi) =o) 6) ee ee (62) 


For n = 3, Eq. (4) is used in addition to Eq. (12) for /=1, 2,3 assuming the user 


altitude A is known. 
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After the position solution is obtained, the control station relays the user position 
to the user via a satellite link as in System I. The frequency assignment and bandwidth 
allocated for each of the satellite links listed below are specified by the Federal Com- 
munications Commission in [Ref. 2]. 


e User-to-Control Station 
User-to-Satellite Uplink 
Frequency 1610 - 1626.5 MHz 
Bandwidth 16.5 MHz 
Satellite-to-Control Station Downlink 
Frequency 5150 - 5166.5 MHz 
Bandwidth 16.5 MHz 


¢ Control Station-to-User 
Control Station-to-Satellite Uplink 
Frequency 6524.5 - 6541 MHz 
Bandwidth 16.5 MHz 
Satellite-to-User Downlink 
Frequency 2483.5 - 2500 MHz 
Bandwidth 16.5 MHz 


¢ Control Station-to-Satellite Command Link 
Frequency 6526.5 - 6541.5 MHz 
Bandwidth 15 MHz 


B. USER-TO-CONTROL STATION LINK 

1. Uplink Bandwidth 

The frequency bandwidth allocated for the user-to-control station data link is 

16.5 MHz from 1610 MHz to 1626.5 MHz. Within this bandwidth, the user must 
transmit a PN ranging code, an identification code, general information, and user alti- 
tude in the three satellite configuration. The user will transmit a standard navigation 
1,023 chip Gold code sequence as the ranging code with the other data. The calculations 
for the user-to-control station are completed as was discussed in System I, utilizing Eqs. 
(5) - (8). System II design parameters are the same as System I. The total carrier-to- 
noise ratio less the 2 dB margin will support data rates of 57 kbps and 18 kbps for re- 
quired signal-to-noise ratios of 15 dB and 20 dB, respectively. The link analysis is 
tabulated in Table 7. 

2. User-to-Control Station Data 
As in System I, the user must transmit a 10 to 17 bit user identification number, 20 bits 
for additional information, and 20 bits for altitude if only three satellites are used. Using 


a standard modulation scheme, this information would be modulated by QPSK and 


BY 


Table 7. SYSTEM II USER-TO-CONTROL STATION LINK ANALYSIS 


Uplink 
Carrier-to- Noise 
Density Ratio 


Downlink 
Carrier-to- Noise 86.6 dB-Hz 
Density Ratio 


Total 
Carrier-to-Noise 64.57 dB-Hz 
Density ne "aed 


Mk a 


If Required ~ === 15 dB 20 dB 


Then Data RIK can be 57 kbps 18 kbps 
































64.6 dB-Hz 








spread with the 1023 chip Gold code. The number of CDMA channels that the system 
will support is the same as in System I. 
3. Satellite-to-Control Station Downlink 
The frequency bandwidth allocated for the satellite-to-control station downlink 
of the user-to-control station data link is 16.5 MHz from 5150 to 5165.5 MHz. The 


satellite functions as a relay only as in System I. 


C. POSITION DETERMINATION 
The position of the user is determined in much the same manner as was done in the 


System I design. The main difference is that there are two distances that must be ac- 
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counted for in the calculation. In System I, the satellites transmit the ranging code and 
the differential distances between the satellite and the user had to be accounted for in 
position error. In System II, the user transmits the code so the user-to-satellite and 
satellite-to-control station differential distances must be taken into consideration when 
calculating the users position. The differential distance measurements are made using 
levels (S805 

The central control station-to-satellite differential distance error does not contribute 
to the user position error as greatly as do the satellite-to-user distance differentials. The 
central control station would be suitably located away from the equator to minimize the 
control station-to-satellite differential distance measurement errors that are inherent to 
the system as was discussed in System I. In the three satellite configuration, where al- 
titude must be entered for position determination, the central control station knows its 
exact altitude, so the only altitude error entering into the position determination is from 
the user. In the four satellite configuration, satellite position errors between the user- 
to-satellite and central control station-to-satellite do increase the user position error but 
not as greatly as first thought. After conducting a computer simulation model run of 
1000 iterations with randomly varying satellite position errors between + 20 and - 20 
meters as was conducted in [Ref. 1]: approximately a 1 to 2 meter additional error was 
added to the overall position error of System I] in companson with System | user posi- 
tion error. It was originally thought that a doubling of the position error would occur 
since an error in the satellite position affects both the user-to-satellite and the central 
control station-to-satellite paths, but was not the case. For this reason, the System I] 


design is a viable alternative to System | design. 


D. CONTROL STATION-TO-USER DATA LINK 
1, Uplink Bandwidth 

The System II control station-to-user data link is designed the same as the data 
link in Svstem ]. The big advantage in System I] is that the bandwidth allocated for the 
satellite-to-user downlink is used solely for the control station to user data link. The 
total bandwidth of 16.5 MHz can be utilized for data transmission because there is no 
ranging code to be transmitted by the satellites. In System I, the ranging code occupied 
10.23 MHz of the 16.5 MHz bandwidth allocated for the satellite to user link, leaving 


only 6.27 MHz for the control station-to-user data link. 


ey) 


2. Control Station-to-User Data 
As in System I, the user ID, latitude, longitude, altitude and general information must 
be transmitted to the user. The System I TDM frame structure will also be utilized in 
System II. In System I], the number of users that can be serviced each second could 
be increased even further due to the 10.23 MHz increase in the available bandwidth. 
The increase in the number of users will depend on the protocols used and is beyond the 
scope of this thesis. The link calculations are completed using Eqs. (5) - (8) and System 
I] link parameters. The values obtained in Table 8 are identical to those in System I for 
the same link parameters. Table 8 is furnished to present the fact that the link calcu- 
lations do not change. 

3. Satellite to User Downlink 

The satellite to user downlink occupies the full bandwidth of 16.5 MHIz from 


1610 to 1626.5 MHz. The satellite operates as a passive relay only. 
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Table 8. SYSTEM Il CONTROL STATION-TO-USER LINK ANALYSIS 


Uphnk 
Carrier-to- Noise 
Density Ratio 


Downlink 
Carrier-to- Noise 
Density Ratio 


Total 
Carrier-to- Noise 
Density Bauo 


ee ee 


If Required <= 15 dB 20 dB 


Then Data icc can be 42.2 kbps 13.3 kbps 






























72.6 dB-Hz 












66.25 dB-Hz 









61.73 dB-Hz 
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IV. SYSTEM COMPARISONS AND CONCLUSIONS 


System I] is the simpler system in many ways. A System II user sends the ranging 
code to the control station so the time of transmission does not require any synchroni- 
zation. The satellite does not require an on-board clock or a code generator and trans- 
ponder to send a ranging code to the user. There exists the possibility of using existing 
satellites or sharing future satellites. The user equipment is also simpler in that the user 
does not have to make time measurements of the satellite-to-user time differentials. By 
reducing the complexity of the user equipment, the cost to each user would also be re- 
duced. 

The number of users that could simultaneously be serviced by System II is also 
greater because the available bandwidth of the control station to user is dedicated totally 
to the link. In System I, the satellite-to-user bandwidth must also support the trans- 
mission of the ranging code from the satellite to the user reducing the bandwidth of the 
contro] station-to-user link. Although there is a slight increase in position error that is 
inherent to System I], the position accuracy is well within expected tolerances required 
of many civilian users. 

This thesis has proposed the concepts of two Global Satellite Location and Tracking 
Systems designed for the civilian user. The advantages of the systems are the low cost, 
user equipment simplicity, and accurate position determination. The user systems are 
cheaper because the computational hardware and software are maintained at the central 
control station. This eliminates the need for installing an expensive computer in each 
user transceiver. 

System I] provides many possibilities for future design due to the fact that satellites 
already in space could be used as communications relays for the system. This eliminates 
the need to design and launch special satellites into space for location and tracking 
purposes only. 

Geostationary Satellite Location and Tracking Systems are viable systems for civil- 


lan position determination. 


Be 
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